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1 Introduction 

Heavy quarks production is an important testing ground for quantum chromodynam- 
ics (QCD), because QCD calculations are expected to be reliable if a hard scale is 
present in the process. In heavy quarks production a hard scale is provided by the 
quark mass. Moreover heavy quarks production can give direct access to the gluon 
density in the proton due to the fact that it proceeds, in QCD, almost exclusively 
via photon-gluon fusion, where a photon from the incoming electron interacts with 
a gluon in the proton giving an heavy quark- ant i- quark pair. Results will be shown 
both for deep-inelastic scattering (DIS), where the virtuality of the exchanged boson 
Q 2 is large, and photo-production, where the Q 2 is equal to zero. Various experi- 
mental techniques are used in order to select charm and beauty events, ranging from 
the measurement of D* cross section to impact parameter analyses. The results are 
found to be compatible with the predictions of perturbative QCD (pQCD). 

This paper is organized as follows. The relevant features of the HERA collider and 
of the HI and ZEUS detectors are described in section [21 In section [31 an introduction 
to the physics of heavy quarks production in ep collisions is given. The sections H] and 
illustrate the tagging methods and the experimental results for the charm quark, 
while E] and [7] do the same for the beauty quark. The charm and beauty structure 
functions are presented in section [HJ The results obtained for the gluon polarization 
by the COMPASS Collaboration are described in section [9j Finally the conclusions 
are drawn in section [TUl 



L On behalf of the HI and ZEUS Collaborations. 
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2 The HERA collider and its two multipurpose 
experiments: HI and ZEUS 

HERA is the first ep collider and consists of two separate rings of circumference 6.3 
km, one a warm magnet electron (or positron) ring with maximum energy 30 GeV and 
the other a superconducting magnet proton ring of maximum energy 920 GeV. The 
rings are brought together at four intersection regions, two of them are occupied by 
the experiments HI and ZEUS. The HERA life can be divided in two parts: HERA-I 
from 1992 to 2000 and HERA-II from 2003 to the middle of 2007. During the first 
period, HERA worked with beam of 27.5 GeV while the energy of the proton beam 
was raised from 820 GeV to 920 GeV. The beam spot had the dimension of 150x30 
(Urn 2 and the integrated luminosity collected by each experiments was about 130 pb _1 . 
At the end of 2000 there was a long shutdown, in which both HERA and the two 
experiments made important upgrades. In 2003 HERA started its functioning and 
it is to foreseen to work up to the middle of 2007. During the HERA-II period the 
energies of the lepton and proton beams remained unchanged: 27.5 GeV and 920 GeV 
respectively. A reduced beam spot (80x20 /im 2 ) and more reliable beams operation 
have enhanced a lot the delivered luminosity to the experiments (~ 180 pb _1 per 
experiment from 2003 to 2005). In Fig. [T]the integrated luminosities per period are 
shown as function of the day of the run. 



HERA delivered 




days of running 

Figure 1: The integrated luminosity delivered by HERA, subdivided into HERA-I 
period and HERA-II, versus the days of running. The HERA-II period is further- 
more divided in HERA-II with electrons (HERA-II e~) and HERA-II with positrons 
(HERA-II e + ). 

The HI [1] and ZEUS j2] detectors are general purpose detectors with nearly her- 
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metic calorimetric coverage. They are designed in order to investigate all aspects of 
high energy ep collisions. In particular both the scattered electron and the hadronic 
system in a hard ep interaction are measured. They are differentiated principally 
by the choices made for the calorimetry. The HI collaboration has stressed elec- 
tron identification and energy resolution, while the ZEUS Collaboration has put its 
emphasis on optimizing calorimetry for the hadronic measurements. The detector 
designs reflect these different emphases. The HI detector has a large diameter mag- 
net encompassing the main liquid argon calorimeter, while the ZEUS detector has 
chosen a uranium-scintillator sampling calorimeter with equal response to electrons 
and hadrons. Cross sectional view of the HI and ZEUS detectors are presented in 
Fig. [2] and [3] respectively. 




Figure 2: Cross sectional view of the HI detector. 



3 Production of Heavy Quarks in ep collisions 

In pQCD, at leading order (LO), two distinct classes of processes contribute to the 
production of heavy quarks (charm and beauty) in ep collisions at HERA. In direct- 
photon processes (Fig. Hk), the photon emitted from the electron enters the hard 
process jg — > QQ directly. In resolved-photon processes (Fig. Hb toHU), the photon 
fluctuates into a hadronic state before the hard interaction and acts as a source of 
partons, one of which takes part in the hard interaction. Resolved photon processes 
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Overview of the ZEUS Detector 
( longitudinal cut ) 




Figure 3: Cross sectional view of the ZEUS detector. 



are expected to contribute significantly in the photo-production regime, in which the 
photon is quasi- real, and to be suppressed towards higher Q 2 . 

Next-to- leading order (NLO) calculations in several schemes are available [3JH]. 
In DIS regime, all approaches assume that Q 2 and heavy quark mass tuq provide a 
hard enough scale to allow the applicability of pQCD and to garantee the validity 
of the factorization theorem. In photo-production regime the hard scale is given by 
the transverse momentum of the heavy quark p ti Q and itlq. In the fixed-order, or 
"massive" , scheme @, u, d, s are the only active flavours in the structure functions 
of the proton and photon. The heavy quarks are assumed to be produced only at 
perturbative level via photon-gluon fusion. This scheme is expected to work well in 
regions where p 2 q ~ roq (if in photo-production regime) or where Q 2 ~ roq (if in 
DIS regime). At higher transverse momenta or Q 2 , calculations based on this scheme 
can break down due to large logarithms ~ l n (Pt,Q/ m Q) (~ ^ n (Q 2 / m Q))- I n this case 
the resummed, or "massless" , scheme B [5] should be applicable. In this scheme, 
charm and beauty are regarded as active flavours (massless partons) in the structure 
functions of the proton and photon and are fragmented from massless partons into 
massive hadrons after the hard process. There are also calculations^ which tempt to 

2 The scheme is often referred to as fixed flavour number scheme (FFNS). 

3 The scheme is often referred as the zero mass variable flavour number scheme (ZMVFNS). 

4 The scheme is commonly referred to as variable flavour number scheme (VFNS). 
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Figure 4: Heavy Quarks production processes in leading order pQCD. 



treat the heavy quarks correctly for all Q 2 . Therefore, at low Q 2 , an heavy quark 
is produced dynamically through the boson-gluon fusion process, whereas, at high 
Q 2 , heavy quark parton densities are introduced. The transition between the two 
extremes is treated in different way by different authors [6]. 

4 Charm production: tagging methods 

The main method used for charm tagging is the identification of the D* mesons using 
the decay channel D* + — ► D q -k^ with the subsequent decay D° — > K~tt + , where ir s 
refers to the low momentum n in the decay. The decay particles of the D* meson are 
reconstructed in the central detector, usually without particle identification. In Fig. [5] 
it is shown a distribution of the mass difference AM = M{Ktttt s ) — M(Ktt) from the 
ZEUS Collaboration. A clear signal is seen around the nominal value M(D*)—M(D°). 
In order to mantain under control the combinatorial background, various cuts are 
made on the p t of the tracks and on the energy of the event. Of course also other 
charmed hadrons were identified and analyzed, such as D + , D s , A c , but with less 
statistics. Finally, the sistematic use of the vertex detectors, first implemented in 
HI and now also in ZEUS, is changing dramatically the perspective of the physical 
analysis in the charm sector as it already happened in the beauty one (see section [6] 
and ED. 
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Figure 5: The distribution of the mass difference, AM = M(Ktt7t s ) — M(Ktt), for 
D* candidates. The D*^ candidates (dots) are shown compared to the wrong charge 
combinations (histogram). The shaded region shows the signal region. The number 
of D* mesons is determined by subtracting the wrong charge background. 

5 Charm production: experimental results 

The status of the charm analysis can be summarized by the two plots of Fig. [61 
where the differential D* cross section as a function of the pseudo-rapidityEI of the D* 
mesons, rj(D*) on the left, and the differential D* cross sections as a function of Q 2 
on the right are shown [7]. 

The plot on the left of Fig. [6] shows the good agreement between the ZEUS and 
HI data. The bands in both plots represent the NLO predictions using the HVQDIS 
program [8], the widths of the bands correspond to the uncertainties in the mass 
of the charm, in the renormalization and factorization scales, in the proton parton 
density functions and in the fragmentation. Rather remarkable is the fact that the 
da/dQ 2 data are well described by NLO calculations over five orders of magnitude. 
Some discrepancies between data and theory are seen in photoproduction: D* pho- 
toproduction cross sections [9] as function of the transverse momentum, pt(D*), and 
tj(D*) show that the predictions from NLO QCD are too low for pt(D*) > 3 GeV and 
T](D*) > 0. Part of this deficit may be due to hadronisation effects. The predictions 

5 The pseudo-rapidity r\ corresponding to a polar angle 9 ( measured respect to the positive z-axis, 
corresponding to the the incoming proton beam direction) is given by rj = -In tan(0/2). 
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Figure 6: On the left, differential D* cross sections as a function of rj(D*), filled points 
from the ZEUS experiment and empty squares from the HI experiment. The bands 
are the NLO predictions of HVQDIS. On the right, differential D* cross sections as a 
function of Q 2 for low Q 2 (dots) and from results on D* production in DIS (triangle) 
compared to the NLO predictions from HVQDIS. The data come from the ZEUS 
Collaboration. 



for single jet and dijet production accompanied by a D* meson should have smaller 
uncertainties from these effects. For that aim the following correlations were studied: 
the difference in the azimuthal angle, A0(D*,jet), between the D* and a jet not con- 
taining the D* meson and those between the two jets of highest transverse energy, 
and the squared transverse momentum of the dijet system, (pr) 2 . For the LO 
2^2 process, the two jets, or the D* and a jet not containing the D* meson, are 
produced back-to-back with A0 = ir and very low px- Large deviations from these 
values may come from higher-order QCD effects. In Fig. [7] the differential cross sec- 
tion as function of the A0(D*,jet) is shown, a large fraction of the produced _D*+jet 
combinations deviates from back-to-back configuration indicating the importance of 
higher order contributions. The available NLO calculations (massive FMNR and 
ZMVFNS [10]) underestimate significantly the observed cross sections in the region 
A(p(D*,}et)< 120°. The cross section da/dAcf)" , see Fig. [HI is reasonable reproduced 
by the NLO predictions in the direct-enriched region, that is x° bs > 0.7d§ , although 
the data exhibit a somewhat harder distribution. In the resolved-enriched region, 
x° bs < 0.75, the data exhibit a harder spectrum than for x° bs > 0.75. The NLO 
prediction of the cross section for x° bs < 0.75 has a significantly softer distribution 
compared to the data. The low- x obs region is more sensitive to higher-order topolo- 
gies not present in the massive NLO prediction. The predictions from PYTHIA MC 

e x° represents the fraction of the photon momentum partecipating to the hard scattering. 
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[IT] reproduce neither the shape nor the normalisation of the data for low and high 
x° bs . However, the predictions from the HERWIG MC [12] give an excellent descrip- 
tion of the shapes of all distributions, although the normalisation is underestimated 
by a factor of 2.5. The fact that a MC programme incorporating parton showers 
can successfully describe the data whereas the NLO QCD prediction cannot indicates 
that the QCD calculation requires higher orders. Matching of parton showers with 
NLO calculations such as in the MC@NLO programme [13J, which is not currently 
available for the processes studied here, should improve the description of the data. 
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Figure 7: _D*+jet cross sections as function of A0(D*,jet) compared with the predic- 
tions of the NLO calculations FMNR and ZMVFNS. 



6 Beauty production: tagging methods 

The HI and ZEUS Collaborations have presented measurements in which the events 
containing beauty are identified in the following manners: using high p? leptons 
(mainly muons) from semileptonic 6-decays, or using the impact parameters of all 
tracks coming from secondary decay vertices (inclusive lifetime tag analysis), or finally 
using double tagged events (D* + //, /i/i). 

In the first method, the transverse momentum of the muon with respect to the 
axis of the associated jets exhibits a much harder spectrum for muons from 6-decays 
than for the other sources. Sometime, in order to enhance the signal to noise ratio 
also the signed impact parameter 5 of the muon track with respect to the primary 
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Figure 8: Cross section for the process ep — > e + D* + jj + X separated into (a,c) 
direct enriched (x° bs > 0.75) and (b,d) resolved enriched (x° bs < 0.75). The data 
(solide dots) are compared (a,b) to the massive QCD prediction with (solid line) and 
without (dotted line) hadronisation corrections applied. The theoretical uncertain- 
ties (hatched band) come from the change in scales simultaneously with the change 
in charm mass. The beauty component is also shown (lower histogram). The data 
are also compared (c,d) with HERWIG (solid line) and PYTHIA (dashed line) MC 
predictions multiplied by the indicated factors. The data come from the ZEUS Col- 
laboration. 



event vertex is used, this quantity reflects the lifetime of the particle from which the 
muon decays. The relative contributions from b, c and light quarks are determined by 
a fit to the p™ 1 distribution or to a combined fit to the p™ 1 and S distributions using 
the shapes of Monte Carlo b, c and light quarks distributions as templates. 

In the second method, the track selection requires full silicon vertex detector in- 
formation. From the measured impact parameter 5 a lifetime significance S = 5/(7$ 
is calculated. Two independent distributions are constructed. S\ is the significance 
distribution of tracks in events with exactly one selected tracks. S2 contains the sig- 
nificances of the tracks with the second highest significance for events with two or 
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more selected tracks. Events in which the tracks with the first and second highest 
absolute significance have different signs are removed from the S2 distribution. The 
subtracted significance distributions are obtained by bin-wise subtraction of the num- 
bers of entries on the negative side from those on the positive side. The subtraction 
method substantially reduces the systematic uncertainties due to track and vertex 
resolutions. The relative contributions from b, c and light quarks are determined 
from a fit to the subtracted S% and S2 distributions and the total number of events, 
using the shapes of Monte Carlo b, c and light quarks distributions as templates. 

In the third method, doubled tagged events, events are selected containing at 
least one reconstructed D* and at least one muon, D* + //, or two muons in the 
final state (/i/i). In order to suppress the various types of backgrounds the charge 
and angle correlations of the D* with respect to the muon and of the two muons are 
exploited. These double tagged measurements extend to significantly lower centre- 
of mass energies of the bb system than measurements based on leptons and/or jets 
with high transverse momentum. Furthermore, these double tagged events permit 
to test higher order QCD effects. For instance, in the photon-gluon rest frame the 
angle between the heavy quarks is 180° at leading order, but at NLO it can differ 
significantly from this value due to hard gluon radiation. 

7 Beauty production: experimental results 

Differential measurements from HI and ZEUS are available for beauty production 
in photoproduction and DIS [14] . [IB] using the lepton+jet(s) tag method. Figure 
[9] shows the differential photoproduction cross sections as a function of the muon 
transverse momentum (on the left) and of the pseudo-rapidity for the process ep — > 
ebbX — > ejjfi. The HI and ZEUS data, which are in reasonable agreement when they 
are compared in the same phase space region (see the da/dr]^ plot on the right side), 
are compared to a NLO calculation in the massive scheme [I]. The NLO calculations 
describe the ZEUS data well. Comparing with the HI data, the NLO calculations 
predict a less steep behaviour for the da/dpt and is lower than the HI data in the 
lower momentum bin by roughly a factor of 2.5; at higher transverse momenta better 
agreement is observed. In DIS (data not shown), the total cross section measurements 
made by the HI and ZEUS Collaborations are somewhat higher than the predictions. 
The observed excess is pronounced at large muon pseudo-rapidities, low values of Q 2 
and muon transverse momentum. 

As said in section El using double tagged events [HI [T7] it is possible to mea- 
sure the b production up to very low pt values. In Fig. [10] the differential cross 
sections as a function of the muon transverse momentum pj, (plot on the left) and 
the muon pseudo-rapidity r] fl (plot on the right), for muons from b decays in dimuon 
events and restricted to the phase space p% > 1.5 GeV and -2.2< 77^ < 2.5 for both 
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Figure 9: Differential cross sections as a function of the muon transverse momentum 
Pt (on the left) and the muon pseudo-rapidity rf (on the right), for muons coming 
from b decays in dijet events. The two sets of data coming from the HI and ZEUS 
experiments were measured in different phase space regions. The full error bars are 
the quadratic sum of the statistical (inner part) and systematic uncertainties. The 
bands represent the NLO predictions convoluted with their uncertaintites obtained 
by varying the 6-quark mass and the renormalization and factorization scales. 



muons are shown. Very good agreement is observed with the PYTHIA+RAPGAP 
[T8] predictions scaled by a factor 1.95 (histogram). Apart from the normalization, 
the leading parton shower approach yields a good description of the corresponding 
physics processes within the entire accessible phase space. The data are also com- 
pared to the absolute NLO prediction in the massive scheme convoluted with the 
hadronization from PYTHIA MC (shaded band). Again, good agreement in shape 
is observed, with a tendency to underestimate the data normalisation. A potential 
trend for increasing data/theory deviations towards low p t and/or high 77, suggested 
by other measurements as said before, is not supported. 

Exploiting the experimental possibilities offered by its microvertex detector, HI 
has measured charm and beauty photoproduction using events with two or more jets 
at high transverse momentum [19J. In this analysis events containing heavy quarks 
are distinguished from light quark events by the long lifetime of c and b flavoured 
hadrons, which lead to the displacements of tracks from the primary vertex (see 
section [6]). This analysis provides the first simultaneous measurement of charm and 
beauty in photoproduction, extending to larger values of transverse jet momentum 
than previous measurements. In Fig. [TT] the measured differential cross sections for 
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Figure 10: Differential cross sections as a function of the muon transverse momentum 
Pj. (on the left) and the muon pseudo-rapidity 1]^ (on the right), for muons from b 
decays in dimuon events. Data come from the ZEUS experiment. The full error bars 
are the quadratic sum of the statistical (inner part) and systematic uncertainties. 
The data are compared to the NLO QCD predictions (shaded band) and to the MC 
predictions (histogram). 



charm (plot on the left) and beauty (plot on the right) as functions of the transverse 
momentum of the leading jet p\ etl are shown. Both charm and beauty data are 
reasonbly well described in shape both by the Monte Carlo simulations (PYTHIA 
and CASCADE^ [H]) and the NLO QCD (FMNR) calculations. For charm, the 
NLO QCD calculation is somewhat lower than the measurement but still in reasonable 
agreement within the theoretical errors, for beauty the disagreement is slightly higher. 
The MC's predict a normalisation which is similar to that of FMNR. The bulk of 
the disagreement between data and NLO calculation, especially for the beauty, is 
observed in the region of small values of x° bs where the prediction lies below the data. 
Restricting the data to x° bs > 0.85, a significant improvement can be obtained: the 
charm cross sections are in good agreement with the NLO QCD calculation both in 
normalisation and shape, the beauty cross sections are also reasonably well described. 

The major part of the results shown in this section were obtained in the photopro- 
duction regime (Q 2 < 1 GeV 2 ), and they differ greatly due to different experimental 
cuts, different tagging-methods. It is difficult to compare each other and also to 
extract a general message from the comparison between data and NLO QCD calcula- 
tions. In order to overcome these difficulties the various measured cross sections were 

7 The CASCADE program implements the kt -factorisation approach instead of the usual collincar 
factorisation approach. In the 73* — > QQ matrix element, which takes the heavy quark mass into 
account, the incoming gluon is treated off mass-shell and can have a finite transverse momentum. 
The calculations are performed at LO, higher order QCD corrections are simulated with initial state 
parton showers. 
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Figure 11: Differential charm (on the left) and beauty (on the right) photoproduction 
cross sections da / dp\ etl for the process ep — > e(cc or bb)X — > ejjX. The full error 
bars are the quadratic sum of the statistical (inner part) and systematic uncertainties. 
The solid lines indicate the prediction from a NLO QCD calculation, corrected for 
hadronisation effects, and the shaded band shows the estimated uncertainty. The 
absolute predictions from PYTHIA (dashed lines) and CASCADE (dotted lines) are 
also shown. 

translated to 6-quark differential cross sections as a function of the quark transverse 
momentum, da(ep — > bX)/dp h Tl in the pseudo-rapidity range \rf\ < 2. In Fig. [12] 
the so extrapolated differential cross sections are shown and compared with the NLO 
QCD (FMNR) calculations (shaded band). The data are in reasonable agreement 
between them, they tend to be somewhat higher than the predictions, the disagree- 
ment is concentrated at low and medium values of p b Tl at high values there is a nice 
agreement. 



8 The charm and beauty structure functions 

The structure functions more frequently studied (F 2 and xF 3 ) are inclusive objects 
and thus contain contributions from both valence and sea quarks. The HI and ZEUS 
detectors have the ability to provide identification of a particular quark flavour open- 
ing so the possibility of studying the contribution of that flavour to F 2 . This is 
particularly important in the case of heavy flavours, as they are likely produced in 
the hard scattering and not in the subsequent hadronisation of the struck parton. In 
other words very precise theoretical predictions can be done as explained in the sec- 
tion[3l Due to the fact that at order a s heavy quark production in DIS occurs through 
boson-gluon fusion process (see Fig. H]), this process involves the gluon density xg 
directly so it gives an experimental handle on this quantity. 

F 2 CC is calculated from the measured charm cross sections as follows: 
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Figure 12: Differential cross section for 6-quark production as a function of the b- 
quark transverse momentum p h T for 6-quark pseudo-rapidity \r] b \ < 2 and for Q 2 < 1 
GeV 2 , 0.2 < y < 0.8. The various points show results from the HI and ZEUS 
Collaborations using different 6-tag methods. The full error bars are the quadratic 
sum of the statistical (inner part) and systematic uncertainties. The dashed line shows 
the NLO QCD prediction with the theoretical uncertainty shown as the shaded band. 
The continuous line shows the k t factorization predictions from CASCADE MC. 



The cross section for cc is calculated from the D* cross section [21] (extrapolated 
to the full phase space) using: 

1 a(ep -> eD*X) 
a(ep^eccX) = - ^ - — } (1) 

where P(c — > D*) is the probability that a charm quark will produce a D* 
meson (about 25%). As said in the sections H] and El the advent of the micro- 
vertex detectors has permitted to distinguish events containing heavy quarks 
from light quark events by the long lifetimes of c and b flavoured hadrons, which 
lead to displacements of tracks from the primary vertex. Furthermore the results 
can be obtained in kinematic regions where there is little extrapolation needed 
to the full phase space and so the model dependent uncertainty due to the 
extrapolation is small. These measurements were done by the HI Collaboration 
(221. 



Finally _F 2 CC is related to ep — > eccX cross-section by: 
d 2 a(cc) 2na 2 



dxdQ 2 Q 4 x 



((l + (l-l/m cc -^ c ), (2) 
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where the small contribution from F£ c is calculated from QCD, while xF% is 
neglected due to the fact that the measurements are made at small Q 2 . 

In Fig. [13] (plot on the left) all the data about F 2 CC are shown as function of x at 
Q 2 values between 2 and 500 GeV 2 . The various data sets, obtained with different 
techniques, are in good agreement between them. The structure function F 2 CC shows 
a rise with decreasing x at constant values of Q 2 . The rise becomes steeper at higher 
Q 2 . The data are compared to calculations using the recent ZEUS NLO fit [23J, 
in which the parton densities in the proton are parameterized by performing fits to 
inclusive DIS measurements from ZEUS and fixed-target experiments. The prediction 
describes the data well for all Q 2 and x except for the lowest Q 2 , where some difference 
is observed. In Fig. [13] (plot on the right) the ratio F^/F^ is shown as function of 
x at fixed values of Q 2 . The charm contribution to F 2 rises from 10% to 30% as 
Q 2 increases and x decreases. The strong rise of _F 2 CC at low values of x is similar to 
that of the gluon density and thus supports the hypothesis that charm production is 
dominated by the boson-gluon fusion mechanism. 

Using the help of the micro- vertex detector it was possible to measure the structure 
function F| b [22] in a similar manner to those depicted for the _F 2 CC . The measurement 
of the b cross section (and so of _F 2 b ) is particularly challenging since b events comprise 
only a small fraction (typically < 5%) of the total cross section. In Fig. [TU the 
measured F| (by the HI Collaboration) is shown as function of Q 2 . The measurement 
shows positive scaling violations which increase with decreasing of x. The data are 
compared with the variable flavour number scheme QCD predictions from MRST [21] 
and CTEQ [25J at NLO and a recent calculation at NNLO [2S]- The predictions are 
found to describe the data reasonably well. The beauty contribution to F 2 , in the 
present kinematic range, increases rapidly with Q 2 from 0.4% at Q 2 = 12 GeV 2 to 
1.5% at Q 2 = 60 GeV 2 . 



9 Polarized gluon distribution 

In this section, results obtained by the COMPASS Collaboration [2Z] on the deter- 
mination of the polarized gluon distribution Ag using the open charm processes in 
polarized deep inelastic scattering are presented. Formally, one may write for the 
spin of the proton: 

X - = ^AE + Ag+ < L z >, (3) 

where AE is the contribution from the quarks and antiquarks, Ag from the gluons 
and the last term the mean contribution of any orbital angular momentum of the 
constituents. While for AE the situation starts to be solid, the challenge remains 
to measure Ag and < L z >. The COMPASS experiment at CERN is a facility for 
spectroscopy and spin physics using hadron and muon beams from the SpS with a 
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Figure 13: On the left plot, the measured F% c at Q 2 values between 2 and 500 GeV 2 as 
a function of x is shown while on the right plot the measured ratio '/ F<i . Data from 
the HI and ZEUS experiments using different charm tagging are shown. The data 
are shown with statistical uncertainties (inner bars) and statistical and systematic 
uncertainties added in quadrature (outer bars). The curves represent the ZEUS NLO 
fit. 

variety of targets and a range of sophisticated detectors for analyzing the final state. 
With a wide range of particle identification devices, the measurement of Ag through 
the boson-gluon fusion production of cc pairs is a primary aim. This first measurement 
was performed by scattering a positive muon polarized beam at 160 GeV on a solid 
polarized target. COMPASS has searched for D° mesons in the decay D° — > K~n + . 
To reduce background the neutral D's was also tagged by requiring them to come 
from the decay D* + — ► D°tt + . In the measurement there is no reconstruction of 
the D° vertex, all the reconstruction is based on the determination of the invariant 
mass and in the identification of the kaon through the RICH detector. The result is 
^9 1 9 = —0.57 ± 0.41 (stat) at a i value of the gluon equal to 0.15 and at a Q 2 = 13 
GeV 2 ; the systematic error is smaller than the statistical one. In Fig. [TSlthe gluon 
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Figure 14: The measured F% b shown as function of Q 2 for various x values. The inner 
error bars show the statistical errors, the outer error bars represent the statistical and 
systematic errors added in quadrature. The prediction of QCD are also shown. 

polarization Ag/g as a function of x at fixed Q 2 is shown. The points represent the 
present LO analyses of hadron helicity asymmetries (mainly from high pt hadrons). 
The result from open charm obtained by COMPASS is also shown (star symbol). It 
is smaller than - but still compatible with - zero. COMPASS performed a NLO fits to 
the spin-dependent structure function g±(x, Q 2 ) world data. Two about equally good 
solutions for Ag(a;, Q 2 ) were found, one with a positive and one with a negative first 
moment AG. 

10 Conclusions 

In the previous pages, part of the results obtained by the HI and ZEUS Collaborations 
in the field of heavy flavours has been summarized. We have seen that their charm 
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fit with AG>0, MS scheme, Q z =3(GeV/c) z 




Figure 15: Gluon polarization Ag/g as a function of x at Q 2 = obtained by NLO 
QCD fits (bands) and from LO analyses of hadron helicity asymmetries (symbols). 

and beauty data are in satisfactory agreement. In photo-production regime, beauty 
and charm data are in general agreement with the NLO predictions, even if beauty 
data are partially slightly higher. Charm production gives a large contribution to the 
inclusive DIS cross section: it was measured with good precision in a large part of 
phase space, NLO QCD calculations describe the data within accuracy. The first F% b 
measurement was also shown. All the presented results come from the HERA-I period, 
much more will come using all the statistics from HERA-II period. In the polarized 
DIS field the new preliminary result on Ag/g from the COMPASS Collaboration using 
open charm was shown. The measurement, considered the most model-independent 
tool to study gluon-polarisation, still suffers from big statistical uncertainties, they 
will be highly reduced using the large amount of data that COMPASS will collect in 
the near future. 
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